Scaling tests of the improved Kogut-Susskind quark action 
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Improved lattice actions for Kogut-Susskind quarks have been shown to improve rotational symmetry 
and flavor symmetry. In this work we find improved scaling behavior of the p and nucleon masses 
expressed in units of a length scale obtained from the static quark potential, and better behavior of 
the Dirac operator in instanton backgrounds. 
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Because of their large computational requirements, full 
QCD lattice simulations profit greatly from the use of 
improved actions, which allow better physics to be ex- 
tracted from simulations at moderate lattice spacings. 
The Kogut-Susskind formulation of lattice fermions is at- 
tractive for full QCD simulations because the thinning of 
the degrees of freedom reduces the computational effort 
and, more importantly, because the residual unbroken 
chiral symmetry prevents additive renormalization of the 
quark mass, and hence eliminates problems with excep- 
tional configurations. 

The lattice artifacts in the Kogut-Susskind formulation 
are order a 2 , unlike the Wilson formulation which has an 
order a artifact, which can be cancelled by the "clover" 
improvement. A third nearest neighbor coupling intro- 
duced by Naik cancels order a 2 violations of rotational 
symmetry in the free quark propagator jl| , and this has 
been shown to lead to an improvement in the rotational 
symmetry of meson propagators 0. Another lattice ar- 
tifact is the breaking of flavor symmetry, signalled by the 
fact that only one of the pions produced from the four 
flavors of quarks has an exactly vanishing mass at zero 
quark mass. This flavor symmetry breaking can be un- 
derstood as a scattering of a quark from one corner of 
the Brilloin zone to another through the exchange of a 
gluon with momentum near ir/a ||. Roughly speaking, 
the cure for this problem consists of introducing a form 
factor for the quark- gluon interaction by smearing out, 
or "fattening" the gauge connection in the quark action. 
In our previous works we have investigated the effects 
of different fat link actions on the flavor symmetry vi- 
olation seen in the pion mass spectrum (see also 
Ref. M.) A theoretically attractive action is the "Asq- 
tad" action studied in Ref ||, which includes the Naik 
correction to improve rotational symmetry, fattening of 
the nearest neighbor coupling to cancel couplings to glu- 
ons with any momentum component equal to 7r/a, and 
a term introduced by Lepage |8) which corrects order a 2 
errors at low momentum introduced by the form factor. 



This "Asqtad" action cancels all tree level order a 2 errors 
1- 

Our previous works tested improvements in physical 
quantities that were directly related to the improvements 
in the action. That is, the Naik term is introduced to im- 
prove rotational symmetry in the quark propagator and 
was seen to improve rotational symmetry in the meson 
propagators. The fat link term was designed to reduce 
flavor symmetry breaking and was seen to reduce the fla- 
vor symmetry breaking in the pion mass spectrum. It 
is perhaps less obvious that other quantities will be im- 
proved, and such tests are the subject of this note. We 
have calculated hadron masses and the static quark po- 
tential at different lattice spacings, and we can compare 
the scaling of these quantities with results with other 
actions. We have also investigated the physics of instan- 
tons on the lattice by computing eigenvalues of the im- 
proved and conventional Kogut-Susskind Dirac operators 
on smooth instanton configurations and on semi-realistic 
"noisy" configurations containing an instanton. 

We test scaling by computing hadron masses using a 
length scale determined from the static quark potential. 
In particular, we will use a variant of the Sommer pa- 
rameter defined by r 2 F(ri) = 1.00, where the commonly 
used r is defined by rlF(r ) = 1.65. (The advantage 
of r% is that it is determined at a shorter length scale, 
where the potential can be determined more accurately, 
and so lattice spacings in simulations with different pa- 
rameters can be matched more accurately using r\.) For 
the quenched potential, this is related to ro and a by 
r i/ r o — 0.725 and r\^fa — 0.85. For quenched simula- 
tions with the single plaquette gauge action, we take the 
lattice spacing from the interpolating formula of Guag- 
nclli, Sommer and Wittig 0: 



log(o/ri) = -1.3589 - 1.7139(/3 - 6) 

+ 0.8155(/3 - 6) 2 - 0.6667(/3 - 6) 3 



(1) 



while for quenched simulations with the Symanzik im- 
proved gauge action we fit a similar formula to the string 
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tension results of Collins et al. pi and our own results 
at P = 10/.g 2 = 8.0 



log(a/ri) = -0.970 - 0.840(/3 - 8) 

+ 0.413(/3 - 8) 2 + 0.304(/3 - 8) 3 



(2) 



with an error of around 1%. We use the Goldstone pion 
mass to adjust the quark mass in the various simulations, 
interpolating to the point where m^ri — 0.778. (This 
somewhat arbitrary value corresponds to am q — 0.02 at 
10/ 5 2 = 8.0.) 
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FIG. 1. Rho masses in units of n. The meaning of the 
symbols is described in the text. 

We calculated the quenched hadron spectrum using the 
improved quark action and Symanzik improved gauge 
action at 10/g 2 — 7.4 and 7.75 using 16 3 x 32 lattices 
archived in our earlier work and on 20 3 lattices at 
10/g 2 — 8.0 generated in our current project, which have 
a lattice spaing of a « 0.14 fm. In Fig. [j] we plot the 
vector meson mass in units of r\ versus the squared lat- 
tice spacing for several combinations of gauge and quark 
actions. The bold squares are from our improved quark 
action, with the Symanzik improved gauge action. The 
octagons are from a simulation with the conventional 
staggered action, and the simple one-plaquette gauge ac- 
tion |l2). The diamonds are the conventional staggered 
quark action, with a Symanzik improved gauge action, 
and the crosses are the Symanzik improved gauge action 
with a quark action containing only the Naik improve- 
ment J|Jll[]. We also show some results for tadpole im- 
proved clover Wilson quark actions. The plusses are from 
the SCRI collaboration |to]], using Symanzik improved 
gauge action, and the fancy squares from the UKQCD 



collaboration p3[| , using the one-plaquette gauge action. 
Fig. U is a similar plot of the nucleon masses in units of 
n . In both of these plots the "Asqtad" action shows scal- 
ing behavior that is dramatically better than the other 
actions tested. 
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FIG. 2. Nucleon masses in units of n. The symbols are 
the same as in Fig. [j]. 

Instantons play an important role in the Euclidean de- 
scription of QCD. They provide the solution to the U(l) 
problem Jl4[ ] , explain chiral symmetry breaking [ ft5[ , and 
at least at a qualitative level reproduce the low mass 
hadron correlators 16 Thus it is important that 

the lattice actions we use to simulate QCD approximate 
the continuum behavior at finite lattice spacing. Here 
we study the topological aspects of the improved Kogut- 
Susskind actions and compare with the standard formu- 
lations. The effect of topology on the spectrum of the 
Dirac matrix has been also studied in |fT8| , and in [fl9|| . 
In fl^l , the microscopic spectral density was computed 
and shown that Kogut-Susskind fcrmions are insensitive 
to topology on coarse lattice spacings (a ~ 0.5 fm). In 
ref. jl9) it was shown that at lattice spacing about 0.07 
fm there is a clear separation between topological and 
non-topological modes. 

We first studied the behavior of the low eigenvalues of 
the Dirac matrix in the background of an instanton by 
computing eigenvalues and eigenvectors of — jp\ using the 
Ritz functional technique pOfl . ip\ is the squared Dirac 
operator restricted to even lattice sites. In an instanton 
background y)\ should have two chiral eigenvectors with 
small eigenvalues, each corresponding to two eigenvec- 
tors of p. We computed the eight lowest eigenvalues of 
— p 2 on smooth instantons. The instantons are put on 
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the lattice by discretizing the continuum formula for the 
instanton gauge field pll, and the resulting gauge 
configuration was smoothed by twenty APE smearing 
sweeps. We calculated eigenvectors on lattices ranging in 
size from 4 4 to 16 4 containing an instanton with radius p 
equal to one fourth the lattice size, p — L/4. Since there 
is no QCD dynamics to define a length scale on these 
smoothed lattices, changing p/a can equally well be con- 
sidered to be varying the size of the instanton at fixed 
lattice spacing or varying the lattice spacing for a fixed 
instanton size. We computed the eight lowest eigenvalues 
of — p 2 for the standard KS action, the Naik action, the 
"Fat7" action( all couplings to gluons with momentum 
7r I a set to zero), the "Fat 7" with the Naik term, and 
finally the Asqtad action (with uq = 1 since tadpole im- 
provement has no meaning on smooth backgrounds) 
As expected we found two small eigenvalues, which ap- 
proach zero as the lattice spacing goes to zero. These two 
eigenvalues were degenerate to the accuracy of the com- 
putation. In Fig. H we plot the small eigenvalue A of p 
for all the actions tested as a function of p/a. For all the 
actions, the small eigenvalues go to zero exponentially as 
a/p decreases. Neither the Naik term nor the fattening 
of the link term affects the eigenvalue significantly. (This 
is not surprising — there is no point in smoothing out 
the link when you are already working on a smooth back- 
ground.) However, the Asqtad action, which differs from 
the Fat7 action in that it contains the small momentum 
correction introduced by Lepage, has smaller eigenvalues 
than the other actions. Wc also measured the chiralitics 
^ = ip^sip of the eigenvectors associated with the small 
eigenvalues. Although \ differs significantly from —1, it 
differs from —1 by an amount proportional to a 2 / p 2 . 

To see if these small eigenvalues persist when the QCD 
dynamics is turned on, we "heated" the p/a — 2 instan- 
ton on the 8 4 lattice. This was done by short quenched 
molecular dynamics trajectories, using the tadpole im- 
proved one loop Symanzik gauge action with (3 = 8.0. 
Since we want to introduce short distance structure with- 
out disturbing the long distance topological structure of 
the initial lattice, we used short trajectories — ten molec- 
ular dynamics steps of size 0.02. We ran for 20 such tra- 
jectories, saving the lattice at the end of each trajectory. 
The resulting sequence of 20 lattices was discarded if any 
of the lattices had topological charge different frome one. 
In total we produced 31 different heating runs in which 
the instanton survived. At the end of the 20th trajec- 
tory the average plaquette was 1.92, similar to the aver- 
age 1.86 of thermalized (3 = 8.0 quenched lattices. For 
comparison, we did 24 runs starting from smooth lattices 
containing no instantons. In Fig. |we plot the 3/4 power 
of the averaged product of the eight smallest eigenvalues 
of p 2 for the Q = 1 heating runs, divided by the same 
quantity from the Q = runs. For all of the quark ac- 
tions, uq was kept at one during these heating runs. This 



quantity, "Detg " , is an approximation to the factor by 
which three flavors of massless dynamical quarks would 
suppress such instanton configurations in a full QCD sim- 
ulation. The abscissa in Fig. ||is 3— (J2av Tr-P A ,„), where 
P^ u is the plaquette. This is related to the amount of 
disorder induced by the heating The vertical line marks 
the average plaquette of a thermalized lattice. As the 
amount of disorder increases, the size of the small eigen- 
values in the Q = 1 runs increases, but the eigenvalues of 
the fat link P rise less than those of the Naik and stan- 
dard KS actions. Roughly speaking, Fig. || shows that at 
a w 0.14 fm three dynamical quarks with fat link actions 
suppress instantons by more than a factor of two relative 
to the conventional action. We have also looked at the 
chirality of the would-be zero modes. Again, as the dis- 
order increases the deviation from —1 also increases, but 
the improved action holds up better. After 20 heating 
trajectories with the p — 2a instanton, the chirality for 
the fat actions is about —0.4 while for the standard KS 
action (and the Naik) it is —0.2. Finally, the degeneracy 
of the two would-be zero modes of — p 2 is lifted as dis- 
order is added. The splitting of these modes is related 
to the flavor symmetry breaking and, as expected from 
our previous studies || , is about a factor of three smaller 
for the fat actions than for the standard Kogut-Susskind 
action. 
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FIG. 3. The small eigenvalue of the Dirac operator as a 
function of instanton size for several actions. The octagons 
are the one link action, the squares are the Naik action, the di- 
amonds are the "Fat 7" action, the bursts are the 'Fat7 Naik", 
and the fancy squares are the Asqtad action. 
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FIG. 4. The approximate instanton suppression factor for 
three quark flavors as a function of the amount of "heating" 
applied to the instanton. The meaning of the symbols is the 
same as in Fig. || 

Lastly, we looked at the eight low eigenvalues of — T/)\ 
for the Asqtad and standard Kogut-Susskind actions on 
424 thermalized 8 4 quenched (3 = 8.0 lattices, which 
have a lattice spacing of about 0.14 fm. At this vol- 
ume we have 1 to 2 topological objects of average size 
p = 2a to 3a per lattice, based on the instanton distri- 
bution measurements in pi]] . Typically, when the topo- 
logical charge is non-zero we find small eigenvalues with 
high absolute value of the chirality (0.2 to 0.3) for the 
Asqtad action. In a chirality versus eigenvalue scatter 
plot the topological eigenvalues clearly separate from the 
non-topological. Similar but weaker separation is also 
observed for the standard Kogut-Susskind action. We 
have also checked the validity of the index theorem. For 
the Kogut-Susskind fermions which represent 4 flavors it 
takes the form 4Q = N_ — N + , where iV_ are the num- 
ber of left-handed modes, N + the number of right-handed 
modes and Q is the topological charge. On the lattice this 
relation does not hold exactly for staggered fermions, but 
the quantities 4 Q and N- — iV+ are strongly correlated. 
For the Asqtad action this correlation is 90%, while for 
the standard Kogut-Susskind is 83%. In both cases the 
topological charge was measured after 200 APE smear- 
ings. 

The code used for generating smooth instanton con- 
figurations and for measuring the topological charge was 
written by Anna Hasenfratz and Tamas Kovacs. Com- 
putations were done on the T3E and PC cluster at 
NERSC, on the T3E and SP2 at SDSC, on the NT clus- 
ter and Origin 2000 at NCSA, and on the Origin 2000 
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